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Abstract: The Abelian Sandpile Model, seen as a deterministic lattice automaton, on two-dimensional 
periodic graphs generates complex regular patterns displaying (fractal) self-similarity. In particular, on a 
variety of lattices and initial conditions, at all sizes, there appears what we call an exact Sierpinski structure'. 
the volume is filled with periodic patterns, glued together along straight lines, with the topology of a tri¬ 
angular Sierpinski gasket. Various lattices (square, hexagonal, kagome,...), initial conditions, and toppling 
rules show Sierpinski structures which are apparently unrelated and involve different mechanisms. As will 
be shown elsewhere, all these structures fall under one roof, and are in fact different projections of a unique 
mechanism pertinent to a family of deterministic surfaces in a 4 -dimensional lattice. This short note gives a 
description of this surface, and of the combinatorics associated to its construction. 
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Introduction. Let A be the lattice in dimension 
4, tensor product of two copies of the triangular lattice, 
A = (ei, 62, 63, 64, 65, 66 | I]l<i<3 = Z]4<i<6 = 0 )z- 

Consider the two-dimensional cell complex containing 
all the vertices and edges of A, and, as (oriented) 
faces, the triangles of the two lattices and the par¬ 
allelograms spanned by pairs (64,64), (62,65) and 
(63,66). We choose the orientation such that the cy¬ 
cles (64,62,63), (-64,-62,-63) and (64,-64,-64,64) 
are upward faces, and similarly with (123) —>■ (456) 
and (14) —> (25), (36). We call a surface a connected 
and simply-connected collection of faces in the cell 
complex above, with all upward faces. 

Embeddings of 64,... ,63 in satisfying the fore- 
mentioned orientation constraints correspond to pro¬ 
jections of the 4-dimensional cell complex on a 2- 
dimensional real space, such that surfaces are mapped 
injectively. An example of such an embedding is 
(64,..., 63) = (w^, w'^,a;°,u;®,w^), where = 

(cos sin ^). In this case, surfaces correspond to 
tilings of regions of the plane, composed only of squares 
and triangles of unit sides, and along directions mul¬ 
tiple of 7r/6. These tilings are called square-triangle 
tilings in the literature. 

Any other projection is topologically equivalent, pro¬ 
vided that 64 -b 62 + 63 = 64 -b 65 -b 66 = 0 and the ori¬ 
entation of the faces is preserved. We call valid such 
a projection. The set of valid projections is an open 
portion of an algebraic projective variety. We call de¬ 
generate projections those on the boundary of this open 
set. Under degenerate projections, the image of some 
faces is a segment or a point. 


A seminal work of de Brujin for Penrose-Ammann 
lozenge tilings |dBr81| has first illustrated the possi¬ 
bility that projections of deterministic surfaces from 
a high-dimensional periodic cell-complex could ex¬ 
plain features of two-dimensional aperiodic incommen¬ 
surable tilings. The square-triangle case discussed here 
shows a similar phenomenon. 

Square-triangle tilings have also distinguished prop¬ 
erties, among which is a relation with Algebraic Ge¬ 
ometry, generalising the well-known connection be¬ 
tween lozenge tilings and Schur functions (see e.g. 
[BP13) ). The algebra of Schur functions has ubiquitous 
three-index structure constants cf called Littlewood- 
Richardson (LR) coefficients |LR34) . When the Young 
diagrams A, p,, D are boxed in a rectangle {d — n) x n 
(as is the case, e.g., when they label cells of the Schu¬ 
bert variety), there exists a relation (Poincare dual¬ 
ity) which acts as complementation at the level of 
diagrams, v v, and the LR coefficients are sym¬ 
metric in all three indices if the upper one is comple¬ 
mented, cf ^ =: As shown by P. Zinn-Justin 

[ZJ09| and Purbhoo [PurOSj . the LR coefficients cor¬ 
respond to the enumerations of square-triangle tilings 
over triangoloids whose three sides are built from A, 
pt and jz, respectively. Two degenerate projections of 
these surfaces reduce to portions of the square and of 
the triangular lattice. As degenerate projections trans¬ 
form some faces into segments or points, the bijective 
correspondence is preserved only if extra integer label¬ 
ings, encoding the disappeared faces, are added to the 
resulting structures. These limiting tilings, together 
with the auxiliary labelings, correspond to the origi¬ 
nal Littlewood-Richardson rule |LR34) in the square 
case, and to the Knutson-Tao (discrete) honeycombs 
[KT991IKTW04) in the triangular case. 


1 
















ASM and square-triangle tilings. The purpose of 
this paper is to illustrate another unsuspected feature 
specific of square-triangle tilings, namely of encoding 
the exact Sierpinski structures that arise in the De¬ 
terministic Abelian Sandpile Model. These structures 
have been identified on various regular two-dimensional 
lattices, under various abelian toppling rules, initial 
conditions and deterministic evolution protocols, and 
square-triangle tilings describe them in a unified way. 

The first occurrences of such structures have been 
presented, by the authors, in [PaolSlICPSlOj . while the 
observation of approximated versions of these struc¬ 
tures (reproduced at a coarse-grained scale, but lo¬ 
cally deformed by some 1-dimensional defects) is much 
older |Ost03) . and has first been made, only on the 
square lattice, for the two most natural determinis¬ 
tic protocols: the evaluation of the identity config¬ 
uration in simple geometries |DM941 ILKG901ILR02] . 
and the relaxation of a large amount of sand put at 
the origin, in the (elsewhere empty) infinite lattice 
|DS(109L[H^IMILPS121[LPST3] . 

The ‘universal role’ of the square-triangle tiling, 
in different ASM realisations, should sound surpris¬ 
ing, as the generic projection gives incommensurable 
parallelogram-triangle tilings and does not live on a 
discrete two-dimensional lattice, as is instead the case 
for the sandpile models we consider. What comes out is 
that, in a remarkable analogy with the mechanism dis¬ 
cussed above for the combinatorics of the Littlewood- 
Richardson rule and Knutson-Tao honeycombs, dif¬ 
ferent lattice ASM realisations occur at different “ra¬ 
tional” points in the set of valid projections (and its 
boundary, of degenerate projections). 

As this short paper is within a series, we do not give 
here an introduction to the Abelian Sandpile Model. 
The interested reader can consult the beautiful re¬ 
view by Deepak Dhar |Dha99] . who first established 
a large part of the theory. For aspects of the model 
more strictly related to the features discussed here, the 
reader can refer to the PhD thesis of one of the authors 
[PaolSj . or the shorter papers |CPS10] and [CPS12] . 
Here we will only concentrate on the aspects concerning 
the surfaces in the square-triangle tiling corresponding 
to the exact Sierpinski structures in the ASM. 

The sandpile configurations are height vectors z = 
{zi}, with variables Zi G N associated to vertices i of 
a graph G = {V,E). There exists a notion of stable 
configuration, and a more restrictive notion of recur¬ 
rent one. Transient is a synonimous for non-recurrent. 
There exists a notion of forbidden sub-configuration 
(FSC), and a stable configuration is recurrent iff it has 
no FSC. More generally, a configuration is recurrent 
over W CV(G) if it has no FSC contained within W, 
thus making recurrency a local notion (like instabil¬ 
ity). Local recurrency and instability are dual notions, 
if we set in the wider frame of multitoppling ASM, as 


first shown in [CPS12] . The toppling matrix A encodes 
the dynamics of the sandpile, and determines a subdi¬ 
vision of into equivalence classes. There exists 

exactly one stable recurrent configuration within each 
class. Unstable configurations z can be relaxed to sta¬ 
ble ones, w = TZz. Stable transient configurations can 
be projected to the unique recurrent representative in 
the class, w = Vz. The operators TZ and V correspond 
to find the fixed point of iterated maps, IZn and Vq, 
corresponding to “rounds” of the procedureu 

A number of structures and operations on square- 
triangle tilings can be introduced, that will reproduce, 
under the various projection procedures, the foremen- 
tioned counterparts in the various ASM realisations. 
We dub all these features of the square-triangle setting 
with the “axiomatic” attribute, as the reason for their 
names emerges only when the projection procedure is 
explicitated. Note that we are not able to reproduce 
all the relevant features of the sandpile model. In par¬ 
ticular, we are not able to reproduce the Oi operators 
(nor their counterparts a| defined in [CPS12j '). The 
main things we are able to reproduce are summarised 
by the following list 

• The notion of (ASM-)equivalence of conhgurations 
is trivialised at the axiomatic level: two tilings are 
equivalent if they have the same boundary. 

• The axiomatic notion of FSC correspond to cycles 
in the tiling satisfying certain local rules. 

• We have an axiomatic notion of Vo, consisting in 
a local deformation along the cycles of maximal 
FSC’s (w.r.t. inclusion). 

• Similarly, we can certify that regions encircled by 
certain cycles will undergo a round of relaxation. 
This gives an axiomatic local notion of unstable 
subconfiguration (USC)ll 

• We have an axiomatic notion of TZq, consisting in 
a local deformation along the cycles of maximal 
use’s (w.r.t. inclusion). 

• We have a recursive description of the Sierpinski 
structures at the axiomatic level. As these struc¬ 
tures in the ASM determine the classification of 
patches and propagators in certain backgrounds 
[Paol3| , this induces a corresponding classiheation 
of axiomatic patches and propagators. 

• A choice of vectors ei,... ,65 G and of “masses” 

{TOi23,ni456,TOi4,TO25, Wae} for the five types of 

tiles, induces a notion of density for the patches. 
This allows to state an axiomatic version of the 
Dhar-Sadhu-Chandra incidence formula, first in¬ 
troduced, for the ASM, in [DSC09) . 

^In TZo one can perform at most one toppling per site, in "Pq 
one adds a single frame identity, and then relax. 

^Corresponding to the waves of topplings |IKP94IIIP98| . 
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Fig. 1 . Bottom left: the Sierpinski gasket of index fc = 3 . 
Bottom right: the structure of the patches, and the role of 
the parameter s. Top: the labeling of the light triangles. 

Sierpinski structures. Let s = {sk ,..., si, sq) be a 
finite string of positive integers, and n(s) = ^Si. A 
Sierpinski structure is labeled by a string s, and n{s) 
is its size. Structures of the same size are equivalent. 

An abstract Sierpinski gasket of index k is defined as 
follows. At index 0, it is just a dark upward triangle. 
At index k + 1, it is obtained from the gasket at index 
k by subdividing all dark upward triangles into three 
dark upward and one light downward triangles, all of 
half the side. Light triangles which are there at index 
k, will remain unchanged at all k' > k. A light triangle 
has index k if it first appeared in a gasket at index 
fc + 1. A gasket of index k has 3^ dark triangles, and 
3^ light triangles of index /i, for 0 < ft. < k. See Fig. 1. 

In the sandpile setting, the triangles of the gasket 
will determine polygonal regions filled with a biperi- 
odic patterns, called patches [Ust03) . Patches may be 
recurrent, transient or marginal, depending on their 
behaviour under the burning test (see [CPSlOj L 

In a Sierpinski structure identified by s, all the dark 
triangles correspond to transient patches, of triangu¬ 
lar shape, with a side of Sk unit tiles. Light triangles 
of index ft correspond to polygonal regions filled with 
recurrent patches. These regions have the aspect of tri- 
angoloids with concave sides, the sides being polygonal 
lines composed of — I segments. The packing of 
unit tiles depends in a certain fixed way on the integer 


ft — ft and the variables Sh> for ft' > ft, and has no extra 
freedom, with an exception: starting at the vertices of 
the triangoloids, we can have a band of a patch with 
marginal density, of width s/i — lH The three bands 
meet at a triangular transient patch. 

A transient patch contains a FSC only if “sufficiently 
large”, namely if it contains at least 7 unit tiles, packed 
in a shape <§§3. Thus, a triangle of side up to 3 units 
Hlled with a transient patch, i.e. the shape , may 
still be part of an overall recurrent configuration. This 
has a consequence on our Sierpinski structures: a struc¬ 
ture with label s is recurrent if and only if 1 < s/i < 3 
for all h < k. These are the structures ultimately ap¬ 
pearing in sandpile protocols. 

Each region of the Sierpinski structure is hlled with 
a periodic pattern. The geometry of every region, in¬ 
cluding the number and location of the unit tiles, is 
determined through a recursive procedure. Also the 
shape of the unit tiles, and their content in terms of el¬ 
ementary squares and triangles, are determined recur¬ 
sively. At this aim it is useful to introduce a labeling of 
the regions of the Sierpinki gasket. We label the dark 
upward triangles with words in the alphabet {a, 6 , c}, 
and the light downward triangles with the same word 
as the dark triangle that originated them. When a tri¬ 
angle of label w is split, the three new triangles, in the 
three directions, have labels wa, wb and wc. We also 
give labels to the three external regions of the triangles, 
as a“^, b~^ and c~^. See Fig. 1. 

A triangle with label w has three larger adjacent light 
triangles, in the three directions, that have labels ct{w), 
(l{w) and 7 (w). These three functions can be defined 
as follows. Let Ou,, Pw and 7 ^, the rightmost position 
along w such that, at its right, there are no more a, b 
or c, respectively; let us call w\{, the truncation of w to 
its first £ letters; let us understand that aa~^ = bb~^ = 
cc~^ = 1. Then a{w) = w\a^a~^, and so on. 

Complex tiles arise from the superposition of more 
elementary ones. Only three tiles are indecomposable, 
and must be given as input. These tiles correspond 
to the three square orientations in our square-triangle 
tilings. The corresponding tilings appear outside the 
triangle, at the three sides. The unit tile of label w 
is composed of the superposition of two copies of the 
tiles of labels a{w), (l{w) and Unless w = 1, 

one of these three words has higher degree than the 
other (say a{'w)). In this case, the six tiles do not 
overlap, with the only exception that the two a(w) 
tiles do overlap exactly on a a(Q;(r(;)) tile. If w = 1, no 
tiles overlap. Each tile has 12 special positions along its 
boundary, which determine the translation vectors of 
the recurrent, transient and marginal tilings involving 
it, and the new tile inheritates its owns positions from 
those of the three subtiles. 


®This corresponds to — 1 parallel type-I propagators, w.r.t. 
the definitions in ICPS10IIPaol3l . 
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Fig. 2. The tile associated to w = cba. The interior orange 
lines describe the decomposition into a(w), /3(w) and 'y(w) 
tiles. The overlap, composed of a a(a(w)) tile and two 
light triangles, is in the middle. The triangles outside the 
tile denote the 12 special positions. Here we have u(P) = 
((61), (2616), (43426162), (34), (5343), (16153435)). 


Dual tiles. Our construction in terms of the vectors 
ei, ..., ee has a number of covariances that allow to 
shorten our description 
Ca-covariance (27r/3 rotations): 

(61,62,63,64,65,65) —?► (62,63,61,65,65,64); 

exchange (123) -f-)- (456) (tt/2 rotations): 

(61,62,63, 64, 65, 65) —>■ (64, 65,65, —61, —62, —63); 

central symmetry (tt rotations): 

(61,62, • • • , 65) —>■ ( — 61, —62, • • • , —65). 

We call polygon a closed curve that is the boundary of 
some square-triangle tiling. A polygon P is determined 
by a cyclic sequence in {1,..., 6,1,..., 6}, where 1, 1 
stand for -|-6i, —ei, and so on. We use the shortcuts 
▲ , T, O and < for the polygons (123), (123), (456) and 
(456), respectively. 

A centrally symmetric polygon P is determined by a 
sequence of the form P = {1112 ■ ■ ■ ik'h'h- ■ ■ *fc), where 
i = i. We use the shortcut {iii 2 ■ • ■ ||) in such a case. 

A polygon P is a dual tile if both the triple of poly¬ 
gons (P, A,T) and the triple (P, I>,<|) (in these pro¬ 
portions) tile periodically the plane. We call a tran¬ 
sient / recurrent hex tiling a tiling of the two forms 
above, respectively. 

The three fundamental parallelogram tiles are dual 
tiles. The dodecagon, (162435||), is another example. 
All the tilings associated to dual tiles, except those de¬ 
riving from the fundamental parallelograms, have the 
topology of a hexagonal tiling: each polygon P is neigh¬ 


bour to other 6 P’s. The fundamental triangles are at 
the 6 triple points, with alternating orientations cycli¬ 
cally along each P. 

To each word w as in the previous section can be 
associated a dual tile P(u>), which is centrally symmet¬ 
ric. The three fundamental parallelograms are (411|) = 
P{a~^) and so on. The dodecagon is (162435||) = P(l). 

A pair of polygons (P, Q) is a dual pair if the sextu- 
plet (P, Q, A, T, >, <) (in these proportions) tiles the 
plane. We call a sq-oc tiling a tiling obtained as above. 

Neglecting triangles (e.g., replacing them with Y- 
shapes), the tiling has the square-octagon topology: 
any Q tile is neighbour to 4 P ones, and any P tile is 
neighbour to 4 P’s and 4 Q’s, alternating^ The funda¬ 
mental triangles are at the triple points of the square- 
octagon topology. Each P and Q tile is adjacent to 
8 and 4 triangles, respectively, alternating dark / light, 
and, within dark and light ones, of opposite orienta¬ 
tions. 

For each w, the pairs of tiles [P{a{w)), P{w)), 
{P{P{w)),P{w)) and (P( 7 (w)),P(w)) are dual pairs. 
For example, the dodecagon and any of the fundamen¬ 
tal parallelograms form a dual pair. 

Exceptionally, and analogously to what happens for 
hex tilings, also all pairs of fundamental parallelograms 
are dual pairs, although with a different topology, and 
with no ordering. 

Each tile P = P{w) appears in two hex tilings, three 
sq-oc tilings as ‘octagon’, and infinitely many sq-oc 
tilings as ‘square’. The union of the positions of triple 
points among all these tilings has cardinality 12. These 
12 special positions break the perimeter of the tile into 
open paths, related by the central symmetry. Thus, a 
list of 6 paths, u{P) = {ui,... ,ue), determines simul¬ 
taneously the perimeter and the special positions, and 
P = (lMi6lt22M34u43u55zi5||). 

The recursive construction, at the level of these 
paths, leads to the formulas (completed by C3- 
covariance) 

{ui6u2)w = (ui6u22u3);g(^) 6 iu^lui6u2).y{w) 

{ui)w = iui)a{w) |a(^6)| > 1/3(16)1 

{U2)w = iu2)l3(w) |a(l6 )| < |/3(i6)| 

(ui)w = (U 2 )iu = 0 a(w) = a~^,/3(w) = 

The geometry of these paths is such that: 

• The sq-oct patches based on a (P, Q) dual pair 
may be adjacent to both recurrent and transient hex 
patches, based both on P and on Q, although with a 
restriction on the direction of the (straight) boundary. 

• The hex transient tiling based on P{w) can be adja¬ 
cent to the hex recurrent tiling based on P{w'), if w' 
is a prefix of w. 

This ultimately leads to the consistency of the con¬ 
struction of the Sierpinski structures (see Fig. 3). 


^This fixes who’s who among P and Q. 
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Fig. 3. Four classes of equivalent configurations: Left: three deterministic configurations, of size n — 6. The two on 
top are stable but transient, and the one on the bottom is recurrent but unstable. Applying V and TZ, respectively, 
we obtain our axiomatic Sierpinski structure, (on the right at s = (1,2,2), thus n{s) = 17). The patch structure is 
highlighted by the orange construction lines, showing the same topology of the Sierpinski gasket in Figure 1. 
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